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Abstract 


We present the general flux tube integrand for MHV and non-MHV amplitudes, in planar 
J\f = 4 SYM theory, up to a group theoretical rational factor. We find that the MHV and 
non-MHV cases only differ by simple form factors which we derive. This information allows 
us to run the operator product expansion program for all sorts of non-MHV amplitudes and 
to test the recently proposed map with the so called charged pentagons transitions. Perfect 
agreement is found, on a large sample of non-MHV amplitudes, with the perturbative data 
available in the literature. 



Contents 


1 Introduction 3 

2 The abelian part 5 

2.1 The dynamical part. 5 

2.2 Charged transitions and form factors. 9 

2.3 Consistency checks. 15 

3 Comparison with data 17 

3.1 NMHV Hexagon . 18 

3.2 NMHV Heptagon. 21 

4 Concluding Remarks 26 

A Pentagon transitions and measures 30 

A.l Summary of transitions. 30 

A. 2 Analytic continuation to small fermions. 33 

A.3 Measures. 34 

A. 4 Zero momentum limit. 35 

B The superconformal charge Q and the flux Goldstone fermion 36 

B. l The zero momentum fermion. 36 

B.2 The commutator of superconformal charges. 37 

C A second example: V\ o P 2 ° P 34 38 


2 













1 Introduction 


The main goal of the pentagon operator product expansion (POPE) program [1] is to pro¬ 
vide an explicit representation of any gluon scattering amplitude in planar Af = 4 Super 
Yang-Mills (SYM) theory at finite coupling. With our current understanding of the POPE 
approach, such an answer will come in the form of an infinite sum, akin to a sort of partition 
function, over all possible excitations of the chromodynamic flux tube of the planar theory. 
More precisely, a scattering amplitude involving n gluons is dual to a polygonal Wilson loop 
with n edges which we decompose into n — 3 successive fluxes/squares. In this picture, the 
transition from one flux to the next is induced by a pentagon operator. Hence, the amplitude 
can be identified with a form factor representation of a correlation function of the n — 4 pen¬ 
tagon operators. The matrix elements of these operators are computed following a somehow 
standard integrable bootstrap [1] . We shall refer to the summand arising in this form factor 
representation, as the POPE integrand. This integrand depends, amongst other things, on 
the n — 5 flux tube states which propagate in each of the n — 5 middle squares 




integrand = integrand(\I> (1) ,... , \I> (ri - 5 )) (1) 

Such states are generically multi-particle states parametrized by a rapidity for each of the 
excitations. For instance, we could have, in the i-tli square, a state with two gluons of 
positive hclicity, a bound state of two gluons of negative hclicity, one scalar and a pair of 
fermions: 1 


= {Fi(ui),F 1 (u 2 ),F_ 2 (u 3 ),fi AB (u 4 ),fi c (u 5 ),fi D (u 6 )} . (2) 

A finite coupling solution for scattering amplitudes in this gauge theory hinges on finding 
explicit expressions for the OPE integrand (1) for any possible multi-particle flux tube state 
such as (2). 

As described in detail below, what renders this seemingly gargantuan task feasible is 
a fortunate factorization of the OPE integrand into considerably simpler building blocks 
which we can analyse separately. It is perhaps worth mentioning from the get-go that this 
factorization is by no means obvious. It stands as another wonderful (but mysterious) = 4 
SYM gift. Were it not for it, one would hardly imagine bootstrapping the multi-particle 
contributions with ease. Indeed, except for the first few particles contributions, almost no 

1 Here A, B,C and D are SU(4) R-charge indices and the indices on F indicate the helicity of the gluonic 
excitation. 
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explicit form factor summands for correlation functions in integrable theories are explicitly 
worked out. A notable exception is the 2d Ising model. The unexpected simplicity we are 
encountering in J\f = 4 SYM theory motivates its portrait as the Ising model of Gauge 
Theories. 

The three building blocks into which the OPE integrand factorizes are dubbed the dy¬ 
namical part, the matrix part and the form factor part, 

integrand = (dynamical part) x (matrix part) x (form factors part) , (3) 

with the latter form factors being non-trivial for non-MHV processes only. Here is what we 
know about these building blocks : 

Dynamical part - This is the part of the OPE integrand that is universally present, which 
applies to all cases, MHV or non-MHV, and which treats all flux tube excitations 
on a same footing, regardless of their quantum numbers / R-charges. It is the most 
dynamical component of the integrand, hence its name, and not surprisingly it exhibits 
the most complicated coupling dependence. Its overall form is however extremely 
simple since it is factorized into a product over elementary pentagon transitions linking 
the various flux tube excitations and a product over square measures and Boltzmann 
weights of each excitation. The geometry of the scattering amplitude, in particular, 
only enters through these Boltzmann weights. All the transitions, measures, energies 
and momentum appearing here are also rather universal. Most of them have already 
been spelled out, see e.g. [2-7], and all of them will be summarized in this paper. 

Matrix part - The matrix part takes care of the 577(4) group theoretical factor of the 
integrand. It can only show up when flux tube excitations with R-indices are present, 
and is otherwise totally absent. (It is also trivial whenever there is only one way 
to distribute the R-indices.) This component of the integrand has the distinguished 
feature of being a coupling independent rational function of the particles’ rapidities, 
with no obvious factorization. Taming this group theoretical factor is an interesting 
algebraic problem on its own [8] but is beyond the scope of this paper. 

Form factors part - Lastly, we have the non-MHV form factors. They are only needed for 
non-MHV amplitudes, which are composed of so-called charged pentagon transitions [2, 
6,9]. Luckily, these form factors are not independent objects. Instead, we can construct 
them from their relation to the bosonic (or MHV) transitions with fermionic excitations 
frozen to zero momentum. Applying this logic, we will obtain their expressions for all 
excitations and transitions. The final result can then be tested against perturbative 
results as well as using self-consistency checks such as parity symmetry. 

The main result of this paper is a complete recipe for writing the two coupling dependent 
factors in (3), that is the complete flux tube integrand up to the matrix part. This can be 
seen as 2/3 of the full POPE program set above and is spelled out in section 2. In section 3 we 
perform several perturbative checks of our POPE elements. The comparisons that we have 
performed test both the form factors and the dictionary proposed in [9] between N fc MHV 
amplitudes and charged pentagon sequences in a rather non-trivial way. 
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2 The abelian part 


In this section we present the expression for the abelian part of the POPE integrand. It 
captures by definition what remains of the full integrand after stripping out the matrix part. 
(In some cases, when there is just no matrix part, the abelian part is of course everything. 
This is the case for instance for states made out of gluons or their bound states, which are 
intrinsically abelian.) As explained in the introduction, the abelian part is composed of the 
dynamical and form factors parts, 

abelian = (dynamical part) x (non-MHV form factors part). (4) 

Conventionally, for MHV, only the dynamical part remains. For non-MHV, the latter re¬ 
mains the same, but form factors should be added to the story. These ones are not really 
independent and can be directly derived from suitable MHV processes, as we shall explain 
in this section. 

2.1 The dynamical part 

We start with the main component. This one captures, in particular, the information about 
the geometry, i.e. the cross ratios <7j,Tj,0j of the polygon, and can be written as [1] 

dynamical part = P( 0 |^ (1) ) e -^ (1) h-Hp(* (1) M^(* (1) )«i 

X P(fh l2) |"^ (1 )) /i(\I>( 2 )) e -- B (^' (2) )i'2-*p(^ ,(2) ) o '2+*m(^'( 2 ))02 

’. ' ' (5) 

X P(T' (2) |T' (3) ) H (tE^ 3 )) e -E^ (3) )T3+ip^))a 3 +irn(^))^ 

x P(^ (4) |^ (3) )... , 

where P(\I/), p(^) and m(T r ) are the energy, momentum and angular momentum of the 
multi-particle state T r . We have n — 5 such states in total, in accordance with the number 
of middle squares in the tessellation, and for each of them we have a corresponding square 
measure /u(fE r ), see [1]. Finally, two consecutive squares with multi-particle states $ and T' 
are connected by means of a pentagon transition P($|T f ) or P(T r |$), where the bar stands 
for the state where all excitations are replaced by their conjugate and their order reversed. 
The fact that each other pentagon appears with such reversed states is a direct consequence 
of the alternating nature of the pentagon tessellation as illustrated in figures 1 and 2. The 
alternating signs multiplying the momenta of the states in consecutive middle squares have 
the same origin. 

The factorization observed above is not a surprise and follows from symmetry considera¬ 
tions of the OPE. In contrast, the simplicity of J\T = 4 SYM theory starts to manifest itself 
as soon as we start exploring the multi-particle nature of the various pieces. What happens 
here is that all the above mentioned blocks factorize further into one- and two-particle blocks! 
To describe this factorization we introduce the notation dp with n — 1,..., N to indicate 
the u-th excitation of the multi-particle state T r . Then, the energy, momentum, angular 
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Figure 1: Every two successive pentagons in the POPE decomposition are flipped with respect 
to each other. Namely, if the cusp of one pentagon is pointing to the right then the next one is 
pointing to the left and so on, ■■■—>• V R —> V L —> V R —> V L 


momentum and measure all factorize into their single particle counterparts as 2 

N 

M*) e~ E ^ T±i ^ +im ^ = J] //(*„) exp [-E^ n )r ± ip(* n )a + im{^ n )<j>] , (6) 

n —1 


where the sign ± = (—1) J+1 multiplying the momenta for states in the j'-tli middle square 
is a simple outcome of the conventions mentioned above, see figure 2. It is convenient to 
use a hatted measure /t to denote collectively the measure and the accompanying Boltzmann 
factor, since these ones always come together. With this notation, the factorization we just 
described would simply read = ][[[ n Most importantly, we observe a similarly 

neat factorization for the pentagon transitions into fundamental 2-particle transitions [1,2, 
4,7,10] 3 

El P($n|^n) 

n,m 

n P(®n\$n’) n P(V m \Vm>) 

n>n' m<m! 

The measures themselves are not independent from the pentagon transitions. On the con¬ 
trary, they can be extracted from the decoupling pole present in 2-particle transitions in¬ 
volving identical in- and out- going particles, 



ResP^i^(u|n) 

V=U 


i 


( 8 ) 


2 In [1] the measure part also included combinatorial factors for identical excitations. These factors can 
instead be associated to the summation over the flux excitation that should be done in a way that avoids 
double counting. 

3 Let us stress again that formula (7) only captures the dynamical part of the transition. In the case 
where $ = 0 and 4/ = 0^, for example, we get from (7) that P<j,<i,(Q\u,v) = 1 /P^^ulv) while in [3] we 
had 0|it, u) = l/(g 2 (u — v + 2 i)(u — v + i)) x 1/P^|^(u|u) which differs by a rational prefactor and by 
the factor 1 /g 2 . The former rational factor is interpreted here as being part of the matrix part and thus 
discarded, while the power otl/g 2 is just absent because of our new normalization of Pm (see appendix A). 
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Figure 2: When inserting the resolution of the identity between each pentagon it is desirable to 
relate one sort of pentagon in figure 1 (say the ones with the cusp to the left) to the other kind 
to render things more uniform. To do so, in previous POPE works, a pentagon with a cusp to the 
left was related to a pentagon with its cusp to the right by a “reflection”, which maps the bottom 
to the bottom and the top to the top. In this paper, instead, we relate a pentagon with its cusp to 
the left to a pentagon with its cusp to the right by a “rotation”, which maps the bottom to the top 
and the top to the bottom. In principle both ought be equivalent. In practice, because of annoying 
minus signs inside some square roots in several pentagon transitions, the second convention is 
considerably more convenient as it avoids several ambiguities that would be present otherwise. In 
this new convention, after rotating the pentagon we relabel the associated rapidities as illustrated 
in this figure. This leads to an alternating sign (—1 ) J multiplying the flux tube space variables as 
written in (5). 


We see that to fully describe the dynamical part all we need are the two-particle pentagon 
transitions between any pair of single particle excitations. Most of them were already written 
down in the literature, see e.g. [2-7]. It was found that, for any pair of excitations {T, $}, 
the pentagon transition takes the rather universal form 

Pm^{u\v) = >*$(«, u)-----, (9) 

S^{u,v) 

where Sq,$(u,v) is the scattering phase (in the symmetric channel) for the excitations T and 
and S*<s,$(u,v) is its mirror counterpart. 1 The functions v) are simple functions of 

the rapidities, involving eventually the Zhukowski variables. In table 1 we present our ansatz 
for these functions for all pairs of excitations. 

In the end, to evaluate the transition we still need to take the square-root of (9) which 
poses some ambiguity on the branch choice. This ambiguity can be partially fixed through 
comparison with data or with help of some reasonable normalization conditions as done 
in appendix A. Explicit expressions for all the transitions and measures are also given in 

4 When f is a gluon or scalar excitation, the mirror S-matrix is given by analytically continuing the 

physical one in the standard way, SW$(it, v) = v). For fermions, the lack of a mirror transformation 

u 4 * * 7 renders the relation between the mirror S-matrix and the physical one less straightforward. It involves a 

so-called anomalous mirror transformation as detailed in [3], see e.g. appendix A.4 therein. 
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P<t, F (u\v) = 

1 , 



1 


{u — V + |) ’ 


:/></> (^1^) 

1 


(u — v)(u — v + i) 

(x + y + - g 2 )(x + y~ - g 2 )(x~y + - g 2 )(x~y~ 

~9 2 ) 

P FF (u\v) = 

g 2 x + x~y + y~ (u — v)(u — v + i) 



(x + y - g 2 )(x~y - g 2 ) 


g\/x + x~y{u — v + |) 

g\Jx + x~y(u — v + |) 


P F ^(u\v) = 

(x+y - g 2 )(x~y - g 2 ) ’ 


J-p F (u\v') = 

g 2 x+x~y+y~ (u — v)(u — v + i) 


(x+y+ - g 2 )(x+y~ - g 2 )(x~y+ - g 2 )(x~y- 

~9 2 ) 

|^) 

(xy - g 2 ) 


y/gxy(u — v)(u — v + i) ’ 


P^(u\v) = 

Vgxy 


(xy - g 2 ) ’ 



( 10 ) 


Table 1: Summary of prefactors for all twist-one squared transitions, with x = x(u),y = x(v) and 
x = ^ (u + — 4 g 2 ) the Zhukowski variable. They agree with those found in the literature up 

to minor redefinitions (see appendix A for details). 


this appendix (altogether with the transitions involving bound states of gluons and small 
fermions [3]). 

In concluding, we recall that there are two main dynamical inputs behind these ansatze. 
The most important (and still mysterious) one is the fundamental relation 

P x \y(u\v) = ±Sxy(u,v)P y \ x (v\u) , (11) 

which comes with a minus sign whenever both X and Y are fermionic. Combined with 
the factorization property (7) it guarantees that general transitions fulfill proper Watson 
and decoupling equations. The other essential constraint is dubbed the mirror axiom which 
states that 

Px\y{u^\v) = Py\x(v\u ) , or more generally Px|*i>(« - 7 |v) = Pq?\ x (v\u) , ( 12 ) 

where —7 denotes the inverse mirror rotation and with v a set of spectator rapidities . 5 
Combined with ( 11 ) the mirror axiom can be used to argue for the ansatz (9) as well as to 

technically, X in this equation is restricted to be a gluonic or a scalar excitation. This is because, as 
explained in [3], the mirror rotation for the fermions is of a more exotic type, mapping the fermions into 
higher-twist excitations. 
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solve for the prefactor T . Finally, it is worth stressing that the solution to such bootstrap 
axioms is by no means unique and comparison with perturbative data is therefore crucial in 
backing up our proposals. This shall be discussed at length below. 

2.2 Charged transitions and form factors 

In the OPE framework, one can cover all helicity amplitudes at once by introducing a super 
pentagon transition [9] 

P = V + x A Va + X A X B V A b + X A X B X C V ABC + X W X D V ABC d , (13) 

where x A is a Grassmann parameter in the fundamental representation of the SU( 4) R- 
symmetry group, and where V Al ...A k , or for short, is the so-called charged transition 
transforming in the k-th antisymmetric product. In the same way that arbitrary MHV 
amplitudes can be described by sequences of bosonic pentagons V, one can generate all 
non-MHV amplitudes by gluing super pentagon transitions together, e.g. 

w = pov o • • • o v y +xlx? x? xf p 12.34 oVo---ov +xixixlxt Pi 23 qp 4 ° • • • ° v + ... 

MHV NMHV NMHV 

(14) 

That this is enough information for recovering the many components of the super Wilson 
loop is not a priori obvious, since the x parameters here are attached to the pentagons in the 
sequence, and not to all possible edges of the loop. However, the missing ‘degrees of freedom’ 
are somewhat superfluous since controlled by supersymmetry, and, as explained in [9], the 
X components are as many supersymmetry independent components as necessary to fix a 
general super amplitude. 

The physics that takes place on the flux tube is essentially the same regardless of whether 
some of the pentagons are charged or not. What can possibly differ is the R-charge flow 
throughout the pentagon evolution. Since the R-charge dependence was factored out into 
the matrix part at the very beginning, one can ask if the abelian part proposed in (7) can 
also be applied as it stands to these charged processes. The answer turns out to be positive 
up to a minor modification: the inclusion of the so-called non-MHV form factors, as sketched 
in (4). The need for these form factors is not a novelty and was previously stressed in [2,7] 
from the study of certain components of the NMHV hexagon. To pave the way to our general 
discussion, let us start by reviewing briefly, on a simple example, why these form factors are 
needed at all, or equivalently why is the dynamical part not enough for describing the abelian 
part of non-MHV amplitudes. 

Consider the ^-component 'P 1234 0 V, or equivalently V o T\ 2 m , of an NMHV hexagon. 
From R-charge conservation the excitations allowed on these transitions are the same as in 
the bosonic MHV case V o V. As such, at twist zero we have the vacuum, at twist one the 
positive and negative helicity gluons, F and F, etc. Despite this similarity, one does not 
expect the transitions, integrands, and full amplitude to be the same for the two processes, 
since e.g. the MHV process treats symmetrically positive and negative helicity gluons while 
the non-MHV one does not. (This is also immediately confirmed at weak coupling by looking 
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at the corresponding amplitudes.) At the level of the POPE integrand for a single gluon, 


VoV = 1 + + J + ■■■, 

^1234 °^ = 1 + J 7 ^fr F (u)f(u) + J f(u) + ... , ( 15 ) 

V 0 V 12 M = 1 + J ^jl F {u)f(u) + J ^fl F (u)f(u) + ... , 

this difference follows from the fact that the gluons are either produced or annihilated in the 

presence of a charged transition pM in the NMHV cases and it results in the gluonic form 
factors / and /. In other words, the dynamical part described before must be completed 
with the knowledge of these form factors, which in the present cases simply read [2] 

/(«) = = 1 //(«), ( 16 ) 

9 

where x ± = x(u ± |) with x(u) — (u + \Ju 2 — 4g 2 )/2 the Zhukowski map of the rapidity w. 
Note that such factor cannot be absorbed in the matrix part both because the matrix part 
is, by definition, independent of the coupling and because in this abelian case there is no 
matrix part. 

Our main proposal is that the same structure persists for generic transitions. Namely, 
the effect of charging a pentagon, ignoring the matrix part, is to dress the abelian part 
by elementary form factors associated to each excitation present on the pentagons. More 
specifically, we propose that the charged version of (7) reads 


r i r(r— 4) 

pM($|^r) = g-g- x 


x n (^r 


X P(§ |«r), 


(17) 


where $ = {$ 1 } is the incoming set of excitations at the bottom square of the pentagon, 
^ = {>Ib} the outgoing set of excitations at the top, and r = 0 , 1 ,..., 4 the amount of 
R-charge carried by the pentagon. The rules of the game are extremely simple. The result 
is factorized and for each excitation, or more precisely for each held creating the excitation 
at the bottom or conjugate held annihilating it at the top, we associate a form factor. The 
form factor can thus be thought of as being attached to the held and represents the net effect 
of charging the transition, as illustrated in hgure 3. 

We immediately verify that the general rule (17) properly reduces to (15) in the case 
of a single gluon. Indeed, as mentioned earlier, the form factor / above accounts for the 
difference between a gluon F produced on top of a charged and an uncharged pentagon, i.e. 

/(«) = P^(0|1.)/P„| F (0|«). (18) 

By convention, or equivalently by applying (7) blindly, P(0|w) = 1, while (17) gives us 

?W(0| U )= 5 °x(M«)) 4 xl, (19) 
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pM(^r|$) 




= (/i*) r x p(v i$) x (h*y 


Figure 3: The net effect of charging a pentagon is to dress each excitation on its edges by a corre¬ 
sponding form factor. This one is attached to the field that creates or annihilates the corresponding 
excitation in the bottom to top evolution picture. The picture above illustrates our conventions, 
with a state T at the bottom being created by a field T and a state $ at the top being annihilated 
by the field <F. The form factors associated to these fields are then {h^) r and (h $) r , with r the 
R-charge of the charged pentagon, or equivalently, the number of times we act with a supercharge 

Q. 


that is 

i 

/0) = ( M «)) 4 & M «) = 

Similarly, we would read that 



( 20 ) 


( 21 ) 


The remaining questions are what are the form factors for the other excitations, why is 
the form factor for a composite state a product of elementary ones, and why is (17) valid at 
all? The quick answers are that, due to their expected simple dependence (see (16)), form 
factors are easily extracted from data analysis, as done in [2], and that their factorized form 
is consistent with all the constraints the charged transitions must fulfill, as shall be discussed 
in section 2.3 below. However, one can do much better than that and actually derive the 
rule (17) directly from the uncharged transitions (7). The important observation [9] is that 
charging a pentagon is the same as acting with a supersymmetry generator on one of its 
edges, as we will now explain. 

To start with, we recall that one can view the pentagon transitions as form factors for a 
pentagon operator acting on the flux tube Hilbert space of states 



( 22 ) 


As explained in [9] , all we need to do in order to add a unit of R-charge to a pentagon V is 
to act with a supersymmetry generator Q on the bottom state or, equivalently, with — Q on 
the top state. Here Q is the unique supercharge that commutes with V and is represented 
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by d x on the super loop (14), see [9]. This leads to the relation : 


Pa(*\&) = (&\Va\*) = (&\VQ a \V) 


Importantly, supersymmetry generators are realized on the flux as zero momentum fermions 
[11]. That is, to act on the state |\E f ) with a supercharge Qa, we add to it a fermion and then 
send the momentum of that fermion to zero. This directly links the charged transitions to 
their un-charged counterparts, allowing us to extract all information about the former from 
the latter. 

Let us show more precisely how this works. Consider for example the transition from the 
vacuum at the bottom to a single fermion ^(u) at the top. Such transition is only possible if 
we equip the transition with R-charge as in the charged transition Pa(0\u) = (^>(u)\Va\0) . To 
obtain the latter transition from an uncharged one we can start with a similar fermion ^(v) 
at the bottom, i.e. from P^(y\u). We now wish to take the limit in which the momentum of 
the fermion '0(u) goes to zero, so that it becomes a supersymmetry generator acting at the 
bottom. In appendix B we show carefully how the zero momentum limit should be taken 
and in particular what is the proportionality factor. We find that 

S |0> = ^ |p) = lim ■ < 24 > 

where, following the notations of [3], the ‘check mark’ on top of the rapidity v, i.e. v, indicates 
that the analytical continuation to the rapidity plane neighbouring the zero momentum point 
(reached at v = oo) has been done (see [3] for further details). Using this prescription, as 
well as the large v behaviours given in appendix A.4, we conclude that 

x = P (^j) * ■ < 25 ) 

where the upper label in P M indicates the amount of R-charge or equivalently the number 
of y’s carried by the pentagon. The fermionic creation transition (25) is obviously of the 
type (17) 6 

^(o|«) = x K( u ) x Po\<,m = x M m ) > ( 26 ) 

with 

h^(u) = (g/x) 1/4 (27) 

the form factor for a single antifermion. As expected it shows a simple dependence on the 
rapidity u of the excitation, once expressed in terms of the Zhukowski variable x = x(u), as 
found earlier for the gluons. 

6 Recall that following (7) we are working in a convention where the un-charged creation transition of any 
excitation X is trivial. With the x-labelling this reads Pg°] Y (0|u) = 1. 
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It is not so much difficult to include more excitations in the top and check the factorization 
of the form factors. Suppose, for instance, that we start with a multi-particle transition 
involving fermions and send the momentum of one of them to zero as prescribed by (24). As 
a result of the multi-particle factorization (7), the multi-particle form factors must factorize 
as well.' 

We can see this at work on simple examples, using the same PM procedure as before. 
For instance, we can at no cost consider the same fermion creation transition with a gluonic 
excitation F a (w ) added on the top. This yields 

= ife 1^#") x > ( 28 ) 

where, again, the upper index indicates that this new transition is taking place on top of 
a pentagon carrying one unit of R-charge. It exactly differs from the chargeless transition 
defined through (7) by the form factors. Indeed, using the factorization of the dynamical 
part (7) together with (25), (27) and (A.4), we verify that 

= 9~*^(u)h F _ a (w)P Q fl Fa {0\u,w ), 

with 

(sign a)/ 4 


in agreement with (20) and (21) for a = —1 and a = +1, respectively. 

Instead of adding matter to the fermion 0 at the top, one can imagine replacing it by a 
pair 00 or a triplet 000 (always with the small fermion 0( v ) at the bottom) and hence access 
to the as-yet-unknown 0 and 0 form factors. A proper analysis would require introducing 
a matrix part, whose main role is to project the pair/triplet to the £77(4) channel with one 
unit of R-charge. However, in both cases, the matrix part plays no role as far as the form 
factors are concerned. 8 One can thus proceed without knowing its explicit form and directly 
relate the 0 and -0 form factors to the large v behaviours of the P^u(v\u) and P^(v\u) 
transitions. Using expressions in appendix A.4 one gets the system of equations 

= g 3/8 Jim U ^y^(^) X P^\<j>iy\u)P^{v\w) = 

19 ** 9 (31) 

^(-i-o=T 4 , 

whose only reasonable solution is 

M w ) = l, h$(u) = ( x /g ) 1/4 

7 To see this in full generality one also needs the matrix part [8]. 

®For completeness, we have indeed that the matrix part is oc l/(w — w + 3i/2) and oc 1/ Y\ i< j{ui — Uj + i) 
for the two cases at hand, i.e. for a state 0(w)0(u>) and 0(ui)0(rt2)0(u3) at the top, respectively. Because 
it shows no dependence at all on the rapidity v of the fermion ip(v) at the bottom, it cannot contribute to 
the form factors. 


(32) 


IR 

1=1 


Ui 


= g 3//g lim 


cusp 


dv 


2 ig dp$ 


x n 


2=1 







(29) 
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Figure 4: Leading twist transition for the hexagon component V 1230 V. 1 . In this example, we assign 
the form factor (/z^(zt)) to the bottom (blue) pentagon and (h^(u)) 1 to the top (red) pentagon. 


Equations (27), (30), and (32) finalize our proposal for the charged transitions (17). A 
couple of consistency checks for it will be given in the next sub-section. One easy test can 
actually be run immediately. It comes from the physical requirement that a pair of conjugate 
excitations should decouple on a charged transition exactly as they do in the un-charged case. 
That is, the square limit (8) must remain the same on a charged transition. 11 Including the 
form factors, this condition translates into 


ResP$i$(u|u) hq, (u)h$(v) 

V=U 


i 


(33) 


which enforces 

h${u)h${u) — 1. (34) 

This relation is easily seen to be satisfied. 

I 11 the end our form factors are all simply given in terms of Zhukowski variables. Putting 
them together, for a polygon with n edges, gives us the full form factors part in (3) as 

n— 5 

form factors part = n 9 8 (hvjr(*)) (hvj/(*+!)) ) (35) 

2=0 


where the index i on top of the matter fields refers to the Pth square, with i = 0 being the 
first one at the very bottom and i = n — 4 the last one at the very top, 9 10 while the same 
index i in r l refers to the Pth pentagon transition between states UP*) and v h'* +1 ), with r, 
units of R-charge. Note in particular that each excitation in a given square is assigned two 
form factors - one for each of the pentagons that overlap on this particular square. A simple 
example is depicted in figure 4. 

9 Put differently, the propagation on the square is diagonal, so it cannot be charged. 

10 We recall that the states in the very bottom (z = 0) and very top (z = n — 4) squares are both vacuum 
with measure and form factors equal to one. 
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As a summary, we can now write the OPE decomposition for a polygon with n edges in 
a rather compact form. Up to the matrix part of course, we have 


P o P o . . . o . 


X = 


„ n—5 

in * 

J 1=0 


j-lqr(i) Q 


n(u~ 4 ) 




n—ri -1 


p R / i (\|/( i )|\|>(*+ 1 b 


(36) 


with A" a choice of x component, with r* Grassmann variables y’s in pentagon i , and with 
pfl/L(^(i)|^(i+i)) = p(^(0|^r(i+i)) Qr p(*j7 (j+1) |*j> w ) fo r i even or odd. 

In what follows we shall perform two sort of checks of our proposal. The first kind of 
checks - with which we will conclude this section - are internal self-consistency checks of the 
POPE proposal. The second sort of checks concern explicit comparison against perturbative 
data and are the main focus of section 3. 

2.3 Consistency checks 

In this section we shall present two consistency checks of the non-MHV form factors presented 
above. Namely, we will first see that they are consistent with parity and finally observe that 
they are compatible with the general axioms for multiparticle transitions. 


Parity 


We have proposed in [9] a simple realization of parity within the POPE approach. In a 
given pentagon with r units of R— charge, the action of parity conjugates its A—charge, 
i.e. r —» 4 — r, which means flipping the chirality of the external particles that are being 
scattered, but also the chirality of the flux tube excitations, namely ift {ft and F -h- F. 
This can be achieved by simply flipping the signs of the angles (ft appearing in the square 
propagation factor of the OPE integrand. 

In more concrete terms, parity symmetry establishes the following relation between the 
OPE components 

P O Po ... O P| x = (P O Po ... O P|_y )|^_0 . (37) 

where X is the complement of the component X , namely 


X 


/ n—4 

n e ^ i=i * ixi x 

i =1 


x->x 


(38) 


For example, we can relate an NMHV component with X = XiXiXiXt f° an 5 MHV with 
its complement X = x\x\xlxl ■ ■ ■ Xn- 4 Xn- 4 Xn- 4 Xn- 4 - For parity to be a symmetry of the 
super Wilson loop, the POPE decomposition should be invariant under r —>■ 4 — r together 
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( 39 ) 


with a flip of the chirality of all the flux tube excitations, namely 


„ n —5 

±n> 

J i =0 


9 




(h\j> (/)) 


ri-n+l 


pfl/L^(i)|^(i+l)x 


» n —5 

$IL 

J *=o 


(rj~ 4 )rj 


9 8 (^4r(i)) 


(4—fi)—(4—r-i+i) 


TiW) I vS"(*+!)> 


P R / L (\|> W |\|> 


J 


up to an unphysical relative normalization. Note that the overall power of g is parity invariant 
by itself. Now, the pentagon transitions and measures change at most by a sign under this 
transformation (namely when fermions are involved) but they do not distinguish the chirality 
of the flux tube excitations otherwise. Finally, using the relation (34) we see that the form 
factor part is also invariant and thus parity is nicely satisfied. 


Bootstrap axioms 

Given the set of elementary transitions (9), the multiparticle ansatz (7) solves minimally for 
the Watson, mirror and decoupling axioms that the multiparticle transitions are conjectured 
to satisfy [1-4]. Any other solution to the same bootstrap equations would differ from (7) 
by a so-called Castillejo-Dalitz-Dyson (CDD) factor. The charged transitions introduced 
in (17) and their associated non-MHV form factors are particular examples of non-minimal 
solutions featuring totally factorized CDD factors. Indeed, as we shall now explain, the 
non-MHV form factors are such as to make the charged transitions (17) a valid solution to 
the multiparticle bootstrap equations. 

As already mentioned, the main feature of the prefactor in (17) is that it is totally 
factorized. As such it is oblivious to the ordering of the excitations at both the bottom 
and the top of the pentagon. It then does not affect the Watson equation encoding the 
transformation property of the transitions upon permutation of excitations in the past or in 
the future. 11 A factorized prefactor is also automatically regular and non-vanishing when 
bottom and top rapidities coincide. The decoupling pole and the decoupling condition for the 
multiparticle transitions are thus preserved, as already stressed around eq. (34). It remains 
to check the mirror axiom (12) which was introduced when bootstrapping the pentagon 
transitions involving gluons or scalars. For this mirror relation to be satisfied by the charged 
transitions, the form factor should obey the relation 

h$(u -7 ) = h$(u). (40) 

This relation is immediately observed for scalars while for gluons it follows from the relation 
x ± (u ~ 7 ) = g 2 /x ± {u). 12 In fact, the validity of both axioms directly follows from the relation 

n Note however that the non-MHV prefactor is not invariant under the exchange of a bottom and a top 
excitation. It cannot therefore be seen as a CDD factor for the elementary transitions (9) without spoiling 
their fundamental relations to the flux tube S-matrix (11). This corrects a misleading statement made in 
earlier versions of this paper. Examples of (factorized) CDD factors for the elementary transitions are given 
at the end of the appendix A. These are redefinitions of the individual pentagon transitions that preserve 
the fundamental and mirror relations, eq. (11) and first eq. in (12), respectively. They also preserve the 
decoupling pole (8), when applicable, but renormalize the square measures. 

12 Beware that mirror transformation takes different form for scalars and gluons [2,12]. 
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between the charged transition and a zero momentum fermion (24) (taken to be one of the 
excitations in (12) for example). 


3 Comparison with data 


We would like now to test our proposal (36) for the POPE integrand against available data 
at weak coupling. Historically, this comparison was absolutely instrumental in unveiling the 
general ansatz for the form factors. 

The data we use is extracted from the package [13], which generates non-MHV amplitudes 
at tree level for any number of particles. The same package also yields one loop amplitudes 
but for the purpose of this paper we restrict our attention to tree-level checks only. 

In the POPE we have essentially five different types of elementary excitations, F, 0,0,0, F, 
and fifteen different pairings of them into transitions. These two numbers are in correspon¬ 
dence with the five independent NMHV hexagons and fifteen independent NMHV heptagons, 
respectively. At leading twist, the hexagons essentially probe the measures and the form fac¬ 
tors. The heptagons on the other hand probe all building blocks of the integrand and, in 
particular, the pentagon transitions involving mixed types of particles. In this section we 
will confront our predictions against both hexagons and heptagons thus probing all these 
ingredients at once. 

At the same time, these checks also provide us with a strong test of the map between 
N fe MHV amplitudes and charged pentagons proposed in [9]. This one maps any N fc MHV 
amplitude, as specified by choosing an //-component of the super loop, to a very precise 
linear combination of the OPE friendly ^-components (14) and vice-versa. It is the latter y- 
components that admit a neat OPE interpretation and thus offer direct access to the various 
pentagon transitions. 

To illustrate this point, let us consider a random //-component, say the NMHV heptagon 
component VO -1 ’ -1 ’ 1 ’ 2 ) multiplying the monomial Cabcd (see figure 8. a for the 

convention of the edge labelling). According to our discussion, we expect it not to have an 
obvious OPE expansion and indeed this is precisely what we find. To see it we extract this 
tree-level component 13 from the package [13] as 


evaluateOsuperComponent [{1, 2}, {3},{},{},{}-, {4}, {}] StreeAmp [7,1] 


and, to make it into a weight free quantity, multiply it by the weights ((—l)i) 2 (l)i (2) 2 with 
(i)j the weight of the twistor Z; t in the j th pentagon (see [9] for more details), which we 
can express in terms of OPE variables using the twistors in [2], 11 We denote this properly 
normalized heptagon component as W*-~ 1,-1,1,2 k In terms of OPE variables we have 




(41) 


13 Actually, the package extracts the ratio function component = yp(*j\M)/yp MHV5 but at tree 

level they are the same. 

14 As is often the case, while this is the correct mathematical procedure, the naive evaluation of yyA -1 ’ -1 ’ 1 . 2 ) 
with the twistors in [2] would yield the very same result since, in this case, the weights simply evaluate to 1. 
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which clearly defies any reasonable OPE interpretation! 15 

Of course, this does not mean that we cannot describe this component within the POPE 
approach. On the contrary, as explained in [9] , once using the inverse map we can express any 
component as a linear component of the nice ^-components which in turn we can describe 
at any loop order within the POPE. In this case, using the expression (21) of [9] we would 
find 


W (-l,-l,l,2) =VuoVo V e -<n-'r 2 -T 2 + l 4± + 


V l2 0? 3 0?4 ( e -2<Tl-n-T2-^-^ _|_ e -<Tl-T2 + _|_ e -<Tl-<72-2T2+ifL + ^. 

^1 j_*02 


*01 *02 


*01 i *02 


P\2 ° Pm ° P 


e _ CT i_ ff2+ ^+^ j + 

.Ml __ ____ _l_Ml 


3 — 2 cri—Ti-n — <Ti-tT2—T 2 + V-+i</>2 cri + ^ L 


+ 


P 123 oVoVi (e- 2<71 ~ T2 -^ + + ^ + + 

P\93 O P 4 O P (g — CTl—n—101 _|_ g — CT1—Tl+201 _|_ e -<72-T2+i<l>2 _|_ g—<Tl-ff 2 —Tl -T2+j^l+i02 _|_ 


—2(Ti—2ri —2(71 


— e 


+ 1 ) + 


P 1234 oV oV [e 


— 2cti—T i + 


^1 


+ e” 


.Mi 


(42) 

It is amusing to see how this precise linear combination of ^-components, each with a nice 
OPE expansion, combines into the component (41) without any obvious OPE picture. Con¬ 
versely, and perhaps less trivially, according to [9] , a generic ^-component is a precise linear 
combination of several r]- components. These are the components that we shall directly 
confront against the integrability inspired predictions in the following subsections. 


3.1 NMHV Hexagon 

Due to R-charge conservation, different NMHV components will support different flux tube 
transitions. The hexagon - made out of two pentagons - is the simplest case where we can 
clearly see this at work. For the components P 1234 0 P and P o P 1234 , encountered before, 
the excitations flowing in the middle square should form an R-charge singlet (for example 
we could have the vacuum, any bound state of gluons, a pair 010234, etc.). For P 12 3 o P 4 
(Pi o P 234 ) we should have excitations with the same total R-charge as for a fermion (anti¬ 
fermion), that is the state should be in the fundamental (antifundamental) representation of 

15 To start with, it simply does not depend on the OPE variables T 2 ,(T 2,</>2 at all, as if only the vacuum 
were propagating in the second square of this heptagon. Even if we were to accept that, other puzzles 
would immediately appear when interpreting the large T\ expansion: the first term looks like a twist zero 
contribution - like the vacuum does - but with a non-trivial 04 and (j>i dependence - contrary to the vacuum. 
There is no natural candidate for what anything like this would be. Also suspicious is the fact that only even 
twists show up once we expand the result out at large 74 . Moreover terms in the near collinear expansion 
shout trouble: their dependence in 04 is so simple that they would not have any sensible Fourier transform 
into momentum space. In short: this component is as weird as it could be from an OPE perspective. 
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Figure 5: (a) OPE friendly edge labelling used in this paper (big black outer numbers) versus 
the more conventional cyclic labelling (small red inner numbers) for the hexagon, (b) The five 
components of the NMHV hexagon and the corresponding excitations at twist one. In the first and 
last component we have written the particle appearing first in perturbation theory. 


£77(4). Finally for P 12 ° P 34 we need the same R-charge as for a scalar, i.e. a state in the 
vector representation of the R-symmetry group. These five components form a basis over 
which we can expand any other component [9]. 

We can write explicitly the OPE integrand for these nice components. At leading twist, 
they are (see figures 4 and 5) 


P 1234 0 V 


P 123 0 P 4 


P \2 ° P 34 


Pi 0 P 234 


P ° P 1234 


1 + f — fi Fl (u) (/iF_i(w )) 4 (h Fl (u))° + f — fiF^u) (h Fl (u )) 4 (hp^u )) 0 + ... , 


o r rjoi 

g~* / 2 ^M«)(M“)) 3 M«) + -- - ’ 


(43) 


9 1 f 7^ Mu) (Mu)) 2 (M«)) 2 + • • • > 


o r rjit 

g~ 3 / ( M «)) 3 + -- -. 


!+ f ^ MM) (h Fl (u))° (h F _ i(u)) 4 + f ^ p, Fl (u) (hp^u )) 0 (h Fl (u )) 4 + .. 


Although in the Erst and last line we have the same allowed excitations, the form factors 
break the symmetry between the positive and negative helicity gluons. In particular, the 
first terms in P 1234 0 P and P o P^a appear at tree level and the last terms are delayed to 
two loops, as confirmed from data. 

As thoroughly described in [3], the contour of the integration C for the fermions is over 
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a two-sheeted Riemann surface and can be conveniently splitted into two contributions: the 
large and small sheet contours. The large sheet contour is performed over the real axis 
with a small positive imaginary part whereas the small sheet contour is a counter-clockwise 
half-moon on the lower complex plane. It is then natural to treat these contributions inde¬ 
pendently as coming from a large fermion 'ipL and a small fermion ^ 5 . I 11 appendix A.2 we 
provide the explicit formulae for the analytic continuation of the pentagon transitions from 
the large to the small sheet. 

Let us now study in detail the component P\ 23 ° P4 at leading twist as an illustration 
of a check against data (see figure 4). The OPE integral splits into two terms corresponding 
to the large and small fermion contributions 


3 f du 3 f du 

V 123 °V4=g-± / — /W M ) ( h $ L (u)) 3 hip L (u) + g~* / — fog (u) {h^ s (u)) 3 h^ s (u ). 

^large ^small 

(44) 

In this case the contour of integration C smaU does not enclose any singularity, resulting in a 
vanishing contribution of the second term at any value of the coupling [3]. We are only left 
with the first term that should be integrated slightly above the real axis. Using the explicit 
expressions for the measure and form factor at leading order in the coupling, we obtain 


P123 o P 4 = e- T+i<t,/2 


du — 


n r 


27T sinh(7ru) 


Jliucr 


+ 0( 9 2 


(45) 


According to the map worked out in [9], this component should relate to a component of the 
super amplitude as follows 

p 123 ° p 4 = ((-i)i ) 3 (4)2 , (46) 

where, as before, the pre-factor (i)j stands for the weight of the twistor Z t in the j th pen¬ 
tagon. In order to extract this component we use the aforementioned package by running 
the following line in Mathematica 


evaluateOsuperCoraponent[{1,2,3},{4},{},{}]OtreeArap [6,1] 

still using the twistors given in [2 ], 16 Upon expanding the outcome at leading order in 
the twist and taking into account the weights, we obtain a perfect match validating our 
conjecture for this particular transition. 

The same type of checks are straightforward to generalize to the rest of the components 
or to one loop level using the same package [13]. These probe the expressions for the form 
factors presented in (35) at weak coupling. We have also verified the correctness of our 
conjectures beyond leading twist when also the pentagon transitions start to play a role. In 
the next section, we probe them more directly using the NMHV heptagon. 

16 Note that one should convert between the OPE friendly and the cyclic labelling of the edges. Figure 5. a 
shows both labellings for the hexagon. 
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3.2 NMHV Heptagon 


The NMHV heptagon is the appropriate laboratory to test the pentagon transitions involving 
all possible pairings of the fundamental excitations. In particular, it is the first polygon where 
the transitions between excitations in different squares arise. All the POPE building blocks 
take now part in and consequently they are all scrutinised. 

The heptagon has fifteen independent components [9] represented in figure 6. Five of 
them, in the top line of that figure, can be constructed in a similar manner to the hexagon 
by charging the outermost pentagons. To generate the remaining ten independent compo¬ 
nents, it is unavoidable charging the middle pentagon. According to [9], charging the middle 
pentagon typically involves rather nontrivial linear combinations of the super amplitude 77- 
components. Since we are also interested in testing this map, in what follows we will focus 
mostly on such examples in which the middle pentagon is charged. 

To leading twist each of the fifteen heptagon components probes a different pentagon 
transition, see figure 6. We verified that their near collinear expansions are indeed in perfect 
agreement with the proposals of the previous section once expanded out to leading order 
in perturbation theory. Interesting as they are, the analysis of these cases follow [2] almost 
verbatim and is therefore not particularly illuminating to present it in detail. Instead, in 
this section we will consider a richer example involving multi-particle states in both middle 
squares. A second example can be found in Appendix C. These examples allow one to 
get a good picture of how generic transitions show up at weak coupling. They also probe 
considerably more structures in a very non-trivial way and allow us to stress the important 
role of the so-called small fermions. The matrix part in the first one is trivially equal to 1 
and in the second case it is simple and has been determined before in [3]. 

The examples where more complicated matrices appear were also tested but we leave 
them for a future publication, where the general construction of the matrix part will be 
presented. All in all, we have tested all possible transitions at tree level up to twists three 
in one square and two in the second one and several twist three (in both squares) transitions 
with simple matrix part (in a total of 429 processes). 

It is instructive to present one such example in detail. We will analyse the V o V 123 0 V 4 
component through the POPE lens. As mentioned above, to make things more interesting 
and nontrivial we will look at some high twist contribution involving several particles and/or 
bound-states. To be precise we shall consider the term proportional to 

e -3n-3i0i x g—3 t2+5/2i<^>2 

as illustration. What flux tube physical processes govern this contribution? To answer this 
question it suffices to list everything that has the right quantum numbers to be allowed to 
flow. In the case at hand we are looking at states with twist 3 both in the first middle 
square and the second square. We are searching for states with helicity —3 in the first square 
and +5/2 in the second square. Finally, we have i+charge considerations. For the sequence 
V o V123 0 v 4 we necessarily have an i?-charge singlet in the first middle square and a state 
in the fundamental (4) representation of 577(4) in the second middle square. All in all, this 
information restricts the matter content enormously. 
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Figure 6: Fifteen independent components for the NMHV heptagon with the corresponding excita¬ 
tions at twist one. The components in the first line involve charging the bottom and top pentagons 
only while the second and third line correspond to the remaining ten components where the middle 
pentagon is also charged. As illustrated here, each such component can be used as a direct probe 
of a pentagon transition. 


22 














































































1/(1) charge -2 — § -1 —\ 0 \ 1 § 2 § 3 


Figure 7: Fundamental excitations together with the effective excitation '^ S F 2 (a.k.a. T>^ip) used in 
the example. 


In the first square, for example, the absolute value of the helicity is maximal, equal to 
the twist of the state. This saturation is achieved for purely gluonic states only, see figure 7. 
There are therefore only three possible states in the first middle square, 

|F_ 1 (w)F_ 1 (v)F_ 1 (k;)) , |F-jMF.jW) , |F_ 3 (w)> . (48) 

The first two are multi-particle states and kick in at higher loop orders. The last one, 
corresponding to a bound-state of three negative helicity gluons, is the only one showing up 
at leading order at weak coupling. 

In the second square things are more interesting. With helicity 5/2 there are only two 
possible states we could envisage: 

\F 1 (u)F 1 (v)^(w )), | ip(w)F 2 (u)) . (49) 

At first we could imagine discarding both since they are both multi-particle states; however, 
when fermions are involved the coupling analysis is more subtle. The point is that fermions 
can be either small or large and in the former case they act as sort of symmetry generators [3]. 
As such each of the states in (49) can be split into two cases depending on whether the fermion 
is small or large. In particular, the second state \F 2 (u)i/j(v)) with the fermion evaluated in 
the small fermion domain can be seen as a supersymmetry generator acting on the excitation 
F 2 (u) thus generating a single effective weak coupling excitation - see figure 7 - and as such 
might show up already at leading order at weak coupling. (At the same time the first state 
\Fi(u)Fi(v)if)(w)) with if} being a small fermion would behave as a two particle state and 
thus show up only at higher orders in perturbation theory.) 

In sum, to match against tree level data it suffices to focus on the process 

vacuum —> F_. 3 (-u) —>■ ijj(w)F 2 (v) —» vacuum, (50) 


23 


which is what we turn to now. In this case the matrix part is trivial. Indeed, the R-charge 
index of the fermion is unambiguously fixed once we pick an R-charge configuration for the 
various y’s. In other words, the matrix part in (3) is equal to one and the full integrand 
is just a product of the dynamical part and the form factor contribution. According to (5) 
and (35) these read 


dynamical part = fi F _ 3 (u) p^w) (i F 2 {v) 


p j\ f 3 H 

\u)Pf. 2 \f 3 {v 

l«) 

P^\F- 2 ( w \ 

H P F 2 w{v\ 

w) 


(51) 


form factors part = (hp 3 (u))° (h F _ 3 (u)) 3 (h^(w)) 3 (h F _ 2 (v)) 3 h^(w) hp 2 (v), (52) 

which we simply multiply together to obtain the POPE prediction 


W F _ 


du f dv 


dw 


3—fh 


— — / — / — (dynamical part) x (form factors part), (53) 


2 tt / 2:/r 


2ir 


^small 


where C sma ii is the complex conjugate version of the half-moon contour mentioned below (43), 
hence running in the upper half u plane, in agreement with the alternating conventions in (5). 
(Equivalently, depending on the ± sign in front of ipa in (5) we use an ±ie prescription for 
integration around zero momentum fermions in the corresponding square.) Plugging all the 
building blocks together, we therefore arrive at 

W f _3^ f 2 = e _3ri_3 ^ 1 x e - 3r2+5/2 ^ 2 x (54) 

f du f dv r dw wT(l + iu ) r(2 - iv)T (1 - m + iv) 2iuai . 2i(v+w)a2 Q( 2) 

J 2 tt ./ 2 vr ./ 2 tt 2 i (u 2 + f) (v 2 + 1) ((u - w ) 2 + 1) } ’ 

» * C Bma n 

and verify that, despite the funny fractional powers of the coupling appearing in the individ¬ 
ual ingredients, the resulting integrand has a regular expansion in g 2 and starts at tree level, 
as expected. (This phenomenon is not accidental and is discussed further in the conclusions 
of the paper.) 

One of the three integrals involves a small fermion ip{w) integrated over its corresponding 
small fermion contour C sma u. An important universal property of small fermions is that they 
can always be straightforwardly integrated out (at any value of the coupling in fact). In 
this tree level example we see that the only singularity inside the half moon encircling the 
upper half plane is the single pole at w — v + i. The fermion integral thus collapses into 
the corresponding residue contribution which freezes w to be attached to the rapidity v in 
a Bethe string like pattern. The interpretation of such strings is that the fermion is acting 
as a symmetry generator on the other excitation in this square, the bound-state of gluons 
F 2 (v). The result of this action is an effective twist 3 weak coupling excitation, see figure 7, 
which is described by the Bethe string. 

In sum, after integrating out the small fermion we end up with the integrations in u 
and v for a single effective particle in each square. The resulting integral can then be 
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Figure 8 : (a) The POPE friendly edge labelling used in this paper (big black outer numbers) versus 
the more conventional cyclic labelling (small red inner numbers) for the heptagon, (b) The NMHV 
heptagon process analysed in this section. 


straightforwardly performed leading to the prediction 


y^F. 3 ^i,F 2 


g—3ri—3r 2 —3i0i+5i02/2 

[e 2ai + l ) 3 (e 2<T2 + l ) 3 (e 2 <T 2 + 2<T i + e 2<72 + e 2ff1 ) 5 


X ( e 13<Tl + 2<T 2 _|_ g gllm+4^ _|_ g gl3o-i+4(T2 


+ 10 e 9ffl+6<72 + 35 e llf7l+6,J2 + 28 e 13rTl+6a2 + 10 e 7(7l + 8ff2 + 60 e 9ai+8(T2 + 105 e 11,Tl+8,r2 
+ 56 e 13ai+8a2 + 5 e 5ai+10a2 + 35 e 7cTl+10<T2 + 105 e 9ai+Wa2 + 130 e llai+10a2 + 55 e 13ai+10a2 

+ e 3ffl+12ff2 + 8 e 5rT1 + 12cT2 + 28 e 7fT1 + 12 °' 2 + 56 e 9 ff i + 12ff 2 + 55 e ii*i+i 2*2 + 2 o e 13 <T i+i 2<T2 ). 

(55) 

This example clearly illustrates the importance of checking the integrability against per¬ 
turbative data. After all, it is clearly a tall order to reproduce any result of the complexity 
of (55). According to the proposal in [9], the amplitude V o V 123 0 V 4 can be extracted from 
standard //-components as 



3 d 
d\3 


W 


(5)3 

((1) 2 (2) 2 (3) 2 ) 3 

(5)3 

((1) 2 (2) 2 (3) 2 ) 3 


((1,2.3,-1)^ + <1.2,3. °>^-) S (56) 

((1,2,3, -l) 3 yvK-b-b-bs) + 2) 3) —1) 2 (1,2,3,0) yt^- 1 “W) 


+ (1,2,3, -1) (1,2,3,0) 2 yy ( - 1 ’° A5) + (1,2,3,0) 3 W (0 ’°’°’ 5) ) 


( 57 ) 


Each of these components can be obtained by running similar code lines in Mathematica 
as in the previous example of the hexagon, using the heptagon twistors in the Appendix A 
of [ 2 ], Once we evaluate the brackets and the weights, we expand the result at large T\ and 
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Figure 9: Table of weak coupling excitations of the flux tube. The fundamental excitations are 
sitting on the boldfaced squares. Additionally one can generate a plethora of effective particles, by 
attaching small fermions if s or ffj s to these fundamental excitations, as they act as supersymmetry 
generators. Besides the excitations represented in this table, one could generate other excitations 
by attaching both ip s and ip s to a fundamental excitation. This gives rise to a third dimension not 
represented here, where excitations involving T>. r ^ are sitting. 


7"2 and pick the term proportional to (47). In this way we obtain a perfect match with the 
expression (55)! 

These are formidable checks of the full POPE construction as they are probing, at the 
same time, the map between charging pentagons and charging edges of [9] as well as the 
(weak coupling expansion of) the various elements of the POPE integrand. We performed 
several other checks of this sort (more than a hundred of them) always obtaining a perfect 
match. For completeness, we present another example in appendix C. We also explored some 
higher loop data but our analysis there was much less thorough. It would be interesting to 
push it much further both in higher twists and higher loops. In particular, it would be nice 
to make contact with the very interesting recently uncovered heptagon bootstrap [14]. 


4 Concluding Remarks 

In this paper we put forward the full coupling dependent part of the POPE integrand. This 
is an important step towards the completion of the OPE program stated in the beginning. 

One of the most amusing features of the POPE which renders it so efficient at weak 
coupling is that the number of particles contributing at a given loop order grows very slowly 
with the loop order. 

To be precise we should recall the notion of weak coupling effective particles and small 
fermions. What happens at weak coupling is that on top of any number of fundamental 
excitations - that is gluons and their bound-states, scalars and fermions - we can add 
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arbitrarily many so-called small fermions ^ or ^ which behave roughly as supersymmetry 
generators acting on the fundamental excitations and morphing them into what one calls 
effective particles , see figure 9. Technically the way this works is that the POPE integrals 
over the small fermions can always be integrated out explicitly leading to string-like structures 
attached to the rapidities of the fundamental excitations. These strings are the mathematical 
depiction of the effective particles. 

To estimate the coupling dependence of each POPE process we simply take the POPE 
integrand (2) for a given matter content and look for the leading power of g of each of the 
factors in the dynamical and form factor part (the matrix part being coupling independent). 
When doing this we should pay special attention to the fermions as small and large fermion 
measures, transitions and form factors scale quite differently at weak coupling. At the end 
of the day, the results we End are remarkably simple. Take the hexagon for instance. We 
readily find that each OPE process scales as 

g M 2 +(K- D! ^y-Te(rt+r b )(8-r t -r b ) , (gg) 

where r& are the number of y’s in the bottom pentagon and r t is the number of y’s of the 
top pentagon. (The fact that r t = ry = 0 is equivalent to rt = ry = 4 is the statement that 
for 6 points MHV and N 2 MHV are trivially related by parity to one another. For NMHV 
we have r t — 4 — r&.) Finally K and M contain the information about the matter content. 
We have 

M — Np+N^-^N^+N^+Np and K — Np-\-—N^, — —N^—Np-‘—N 1 j, s -\-—N^ s . (59) 

Here Np (Np) indicate the total number of gluons of positive (negative) helicity or their 
bound-states and denotes the number of usual (i.e. large) fermions while N^ s is the 
number of small fermions. The exponent in (58) is always an even number (which is nothing 
but twice the number of loops). Let us expand a bit on the physics of (58) and (59). 

We see that M is the total number of excitations except for small fermions. In other 
words, since the small fermions can always be integrated out to simply change the flavour of 
the other excitations, M is nothing but the total number of effective excitations, 

M = number of effective excitations. (60) 

Take for example a process with 6 effective excitations for the MHV hexagon. From (58) we 
see that it will first show up at eighteen loops (if K = 0) or even later (if K 0). Up to 
seven loops, for instance, we shall never need more than 3 effective particles to describe the 
six point MHV scattering amplitude! 

Next we have K which contrary to M can be positive or negative. Consider first a 
configuration containing only fundamental excitations and no small fermions. Then K is the 
average of ratios of the 17(1) charge to bare twist of each excitation. It is a sort of measure 
of helicity violation. Equivalently, if we associate a weight of 1 — j/2 to an excitation in the 
j'-th row in figure 9 then K is the average of these weights which we can therefore depict as 
a vertical centre of mass position of sorts in this figure. This second description is the most 
convenient one to generalize to the case where we add small fermions to the mix. As indicated 
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by the red arrows in figure 9, such small fermions transform the various excitations moving 
them up or down in this table. With the precise signs in (59) we see that the interpretation 
as the average of rows is perfectly kept when we move to a description purely in terms of 
effective excitations. In short, 

K = vertical centre of mass of figure 9. (61) 


In sum, OPE processes are coupling suppressed by two effects: large number of particles 
and large helicity violations. 

We see that for a given number of effective particles M the processes that minimize the 
loop order are those for which K = rb ~ n . For the NMHV hexagon this means K = which 
translates into the intuitive statement that the centre of mass position K should coincide 
with that of the leading order excitation flowing in that hexagon as identified in figure 5b. 
Conversely, for a given number of effective particles M the processes that maximize the loop 
order are those for which K is as large as possible (positive or negative depending on weather 
rfe is less or greater than ?y). For MHV, for instance, the processes that maximize \K | are 
the states containing nothing but gluons. In that case (58) reduces to 

g 2Nl+2N% 



reproducing the counting in [4], see equation (44) there. As explained in [4], because these 
states maximize U( 1) charge they are very easily identified in perturbation theory by taking 
a so-called double scaling limit. This sector also appears to be instrumental in making 
contact with the hexagon bootstrap program [15-19] as thoroughly explored recently in [20], 
using technology from [21]. In this regard, it would be interesting to see if the OPE loop 
thresholds (62) parallel some sort of complexity jumps or qualitative changes in the hexagon 
function representations. 

A similar loop counting can be straightforwardly performed for higher n -gons where we 
find that a given OPE sequence first shows up at g 21 with the number of loops l given by 




Ml - M,M i+l + 



- K t K l+ 1 



i =1 


Here Mi and Ki are associated to the i-th middle square and are defined, for a given matter 
content flowing in that square, exactly as above (59 ). 17 Obviously, the interpretation as (60) 
and (61) continues to hold. Similarly, the r* are the number of y’s in the i-th pentagon 
with ro = r\ for i — 1 and r n _3 = r n _4 for the boundary cases. Note in particular that for 
heptagons and highers the loop counting is not simply a sum of squares. Now some terms 
contribute with a minus sign. An interesting outcome of this fact is that with higher n-gons 
we can often engineer processes with large number of particles at relatively low loop order 
by considering polygons with many edges and slowly injecting more and more particles and 
x’s as we move along the tessellation, in the same way as one efficiently starts a car by 
gently pushing it down a road. The longer the road, and the more help one gets along the 
way, the easiest it is to start the car. This is often quite useful if one wishes to test higher 

17 The values for M „_4 and A '„-4 should be set to zero. 








particle contributions at weak coupling without going to prohibitively high loop orders. Some 
examples of this enhancement were already discussed in a purely gluonic context in [4], For 
instance, we can set all r’s to two, start with a scalar in the first square and add one scalar 
in each square until a maximum halfway through the tessellation where we start decreasing 
one by one the number of scalars. A simple counting exercise shows that this configuration 
will first appear at tree-level for any polygon. To give an example, we could check a process 
involving 4998 scalars in a myriagon already at tree-level! 

In [8] the full matrix part for any n-gon, the last missing piece in the POPE integrand, 
is analysed and in [22] the outcome of this analysis shall be unveiled for the hexagon case. 
Together with the result herein presented these works flash out the initial proposal in [1] 
to completion. Nonetheless, in our view, they are not simply the end. Having a fully non- 
perturbative proposal for scattering amplitudes in a four dimensional gauge theory is very 
exciting but more exciting still is the perspective of using it to extract sharp physics out 
of it at weak, strong and finite coupling, thus substantially enlarging our understanding of 
non-perturbative quantum held theory. 
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A Pentagon transitions and measures 


In this appendix we summarize our knowledge about elementary transitions Px\y, with 
(X, Y) being any pair of flux tube excitations. Their general structure is 


Px\y(u\v) = fx(u)F X Y{u,v)f Y (-v) 


x exp 


2 (k x (u) 


ikx(u)Y ■ M ■ k y (v) + 2 i(K X {u) + in x (u)Y ■ M ■ n Y {y) , 


(64) 

where all the objects in the exponent were explicitly given in the appendix C of [3] for all 
sorts of excitations of the flux tube (with J\4 — Q ■ M in the notations of [3], see also [4]). 
We also found convenient to strip out the factor 


log f x (u) 


I j(M2gt) - 1 ) 


\Jo{2gt) + \ - e ~ qxt e~ iut 
e* — 1 


o 


(65) 


for each excitation, with J 0 (z) = 1 + 0(z 2 ) the Bessel function of the first kind and qx = 
—1/2, 0,1/2,1 ,... for scalar, large fermion, elementary gluon, bound state of two gluons, etc. 
(Note that in the case of a small fermion, i.e. X = ij} $ or t/s, we have fx(u) = 1 identically, 
see A.2 below.) 

The factor (65) as well as the term in the exponent above are quite universal and, in par¬ 
ticular, only depend on the absolute values of the (7(1) charges (e.g. they cannot distinguish 
between (X,Y) = (^>, ^>), (-0,^), (0,0), or (0,0))- As suc h) the function Fxy(u,v ) has the 
same conjugation property as its parent transition. Since all our transitions obey 18 


Px\y(u\v)* = P Y] x(v\u ), ( 66 ) 

upon complex conjugation (for real rapidities), with of course 0 = 0 for a scalar, F a = TL a 
for a gluon, etc., then the exact same relation holds true for the corresponding functions F. 

It is also interesting to note that both log fx and the exponent in (64) are of order 0(g 2 ) 
for small g. (This estimate is not uniform in the rapidities and holds only away from the 
locations of singularities, which are at imaginary half integer values at weak coupling.) The 
leading order weak coupling results can thus be directly obtained from the prefactors F XY - 


A.l Summary of transitions 

Knowing the transitions is equivalent to knowing the prefactors F in (64). For them, which 
as we just said are also all we need to know to leading order at weak coupling, we have 
the following lists. (Up to few exceptions involving gluonic bound states, all the transitions 
given below already appeared in [2-7]. 19 At the end of this subsection, we comment on the 

18 We failed to find a reason for this simple property, but noticed that it is consistent both with the 
fundamental relation (11), since Sxy(u,v)* = S Y x(v,u ) = Syx(v,u), and with the mirror equation (12). 
In the latter case, one needs to use that u ~ 7 turns into u +1 upon conjugation and that (12) is equivalent to 
P X \y{u\v 7 ) = P Y \x(v\u). 

19 Mixed transitions involving gluonic bound states were independently obtained by A. Belitsky [23]. 
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‘difference of normalizations’ between our transitions here and those obtained in these series 
of papers.) 


Transitions involving a gluon or a bound state of gluons 


We start by the cases involving a gluon or a bound state of gluons. 

The purely gluonic transitions are given by 

F Fa F b (u,v) = \J (T[ +a V -b ] — g 2 ){x\~ a \y[+ b \ — g 2 )(x[ +a V +b l — g 2 )(xF a ]yl~ b ] — g- 


F F a F b (u,v ) = 


(-l)&r(Md4. + iu - iu)r(^l±^ -iu + iv ) 
c/ 2 T(l + H + iu) r(i + JM - iu)r(l + -iu + iv) 

1 

a;[+ a ]j/[+*r) 

T(1 + + iu - iu) 


, for ab > 0 , 


r(l + ^ + m)r(l + f -iv) 


, for ab < 0 , 


and the mixed ones by 


y , a; [+a] a; [-a]r(i -f M -)- iu — ill) 

Ff ‘ M = = s r(i + M + i U )r(i-») ' 


(67) 


( 68 ) 


for any a, and by 


Ff^{u,v) = Fj, Fa (-v, -u) 


iy(x [+a] x [ a ]) 3/4 T(| +iu- iv)yj(l - ^^)(1 
g 3 / 2 T(l + § + m)T(l - iv) 



) 


J 


F F- a ^(u,v) - F i>F _ a (- 


-u = 


(; 0 +a W a]jl/ 4 p^ + | -|-iu — i V ) 


9 1/2 r(i +1 + iii)r(i - ’ 


(69) 

for a > 0. In all cases, we have x = x(u), = x{u ± ia/ 2), y = x(v), y^ = x(v ± 

ib/2),x(u) = \{u+ \/u 2 — (2 g) 2 ) and T(^) the Euler Gamma function. Transitions involving 
■0 can be obtained by conjugating those with -0, as in (66). 


(Note that all the factors above are normalized such that the associated transitions are 
equal to 1 to leading order at strong coupling, in the perturbative regime, i.e. for excitations 
with momenta of the same order as their masses. This is the expected decoupling property 
of the gluons at strong coupling.) 
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Transitions involving only scalars or fermions 


We proceed with the remaining set of functions Fxy involving scalars and fermions. They 
read 

T (iu — iv) 

HU,V gT(\ + iu)T(\-iv )’ 


F h {u,v) = F^(u,v) = 


y/yT (l + iu-iv) 


gT{\ + iu) T(1 — iv) ’ 


F^(u,v) = —F^(u, v) = 
F^(u,v) = -F^(u,v) = 


(xy) 3 P T(1 + iu — iv) 


g 3 / 4 V{l + 7u)T(l — iv) \J xy — g 2 
i(au/) 1 // 4 r(ra — iv)^/xy — g 2 


(70) 


g 5 / 4 T(l + ra)T(l — iv) 

(The branch choice for the mixed transitions above was mostly driven by the goal of getting 
the sign-free large (positive) v behaviours (83).) 

The lists (67), (68), (69) and (70) cover all the pentagon transitions of the OPE program. 


Comparison with the literature 


As alluded to before, the pentagon transitions listed before already appeared in the litera¬ 
ture. We found useful however to redefine some of them, still preserving the fundamental 
relation (11) and the mirror axiom (12). For instance, comparing the scalar transition in (70) 
with the one appearing in [2], we find that 


F(f,\(f>(u\v) here 9 P<j>\<j>(.'U‘\v) J2] • 


(71) 


Clearly, such an innocent rescaling cannot alter the validity of (11) and (12). A slightly 
more involved redefinition is found when comparing our mixed transitions in (69) with those 
found in [6]. We get 


-fV|y( M KOhere — 


FF\ipi U \ V )here ~ 


Vx + x~ 

9 

9 

V x + x~ 


x P F ^(u\v ) 2 [6] , 
x Pf\4>( u \ v ) 2 [6] ■ 


(72) 


Nonetheless, here as well, both expressions fulfill the same axioms (11) and (12), and thus 
differ only by CDD factors. Along a somewhat similar vein, 20 we have 


,3/4 


(tt|u)here ^ (xy) 4 i 4 ^ Fip\i/){u\v) [3] , P ^|^( 


^PVhere 


^PxP^iiu Mm, (73) 


and also 

P'ijj'ipiF 1 |^)here Pi^ijjiu |^)here • Punjj ( ^ I F) here Ppy ( ^ I here * (74) 

20 In case of fermions, it is much harder to check which choice is better as far as the mirror property is 
concerned, since the mirror algebra is anomalous for fermions and thus harder to implement, see appendices 
of [2] for more details. 
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while the convention of [3] was 

p w( u \ v ) [3 ] = P^(u\v) [ 3 ] , Pj$(u\v) [ 3 ] = P^(u\v) [3] . (75) 

These are all the differences we could find. As we can see, they almost all relate to the 
choice of CDD factors, which (by definition) is the sole freedom left after imposing the 
bootstrap axioms. 21 This ambiguity is not without connection to the freedom we have to 
reshuffle factors in the definitions of the various elements of the POPE integrand, provided 
we have the same physical integrand once we multiply them together. This is precisely 
what happens when relating the integrands that we can produce here with those derived 
in [2,3,5-7] using slightly different transitions, measures, and, sometimes, rules for assigning 
form factors. Our conventions here have the slight advantage that no form factors are 
needed for MHV amplitudes; equivalently, they are normalized such that the relation between 
charged pentagons and zero momentum fermions is most straightforward. It is nonetheless 
clearly desirable to find a better handle on the ambiguities mentioned above. Additional 
constraints should, for instance, follow from the nonlinear nature of the anomalous mirror 
rotation for fermions and/or from certain fusion properties the various transitions could obey. 
It remains to work them out and see how much they can tell. 


A.2 Analytic continuation to small fermions 


It is also convenient to store the representation for small fermions. This one is obtained by 
direct analytical continuation, x(v) —> x(v) = g 2 /x(v), and it was thoroughly exemplified 
in [3]. The squared transitions, for instance, read as : 


P<t>4,{u\v) 2 


P<jjip(^i 7 ) 

(u-v + | )S i "w(u,v) ’ 


Pftp(u\v ) 2 

P F 4 ,(u\v ) 2 

P^nj; (u|7 ) 
P^ii) (^ 1^0 
P^iulv) 2 


g\fx~ 

v x~y(u - v - %)S Fl /,(u, 

v) 

(x+y 

- g 2 )(x~y - g 2 )S^(u 


(x + y 

- g 2 )(x~y - g 2 )S F j(u, 

v) 

g\/ x^ 

'-s- 

•«|(M 

1 

1 

B 

v) 


a/ 9%V Sip4,(u, h) 


(xy - 

g 2 )(u -v + i)S^(u,v 

)’ 

(xy 

-9 2 )S^(u,v) 


yjgxy (u - v)S^(u,v) ’ 



(xy - g 2 )S^(u,v) 



y/gxy(u -v)(u-v + i)S^ 4 ,(u,v) ’ 


y/gxyS^{u,v) 

(xy - g 2 )S^(u,v) ’ 


(76) 


21 In the conventions of [3] the fundamental axiom (11) for was not including a minus sign in front of 
the S-matrix; hence differences in transitions related to this sign are not strictly speaking of CDD type. 
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where, again, ‘check marked’ rapidities indicate analytical continuation to the small momen¬ 
tum sheet. 

More explicitly, and including bound states as well, we can take all the transitions listed 
before and perform the continuation. The explicit form of the analytically continued tran¬ 
sitions still preserves the structure (64), but the prefactors as well as the functions in the 
exponent are changed. The continuation of the exponent was explicitly worked out in [3] 
and here we provide once again only the expressions for the prefactors. We stress in par¬ 
ticular that f^ s {u) 7 ^ /y('h) and F x ^ s (u,v) ^ F x ^(u,v), since upon continuation of the full 
transition some extra terms are produced by the exponent in (64) and transferred to the 
prefactors. Instead, the correct analytic continuation produces /y s (n) = 1 and 


F F a ^ s (u,v) = F^ sFa (-v, -u) = 

Ff^ s (u,v) = F^ sK (-v, -u ) = 
and 

F Hs (u,v) = 1/y/y, 


y/g(u — v + iaf 2 ) 


?/ ( z[ + “W-“] ) i/y (1 - ’ 

' x [+°] x [— °l)i/4 


for a > 0 , 


9 2 \ 9 2 


V9 


for «<o, 


x [+ a \y' 


Fi> s <i>(u, v) = 1 /y/x, 


9 


1/4 


g- 


Fi>ip s (u,v) x l/4y3/4^ \j xv ’ F^ s (u,v) 


1 (X 


xy 


g 1,A \y 


1/4 


g l,i 

F^si>{u,v) = - x in yX iJ \j 1 


g_ 

xy 


i Fibc/ib f) 


1 /y\ 1/4 


Fil>si>s( u i v ) — 


(X‘ 


y ) 1 ' 4 


(/ 1 / 4 (u — V ) 


l- 9 -, 

xy 


i/j s ip s 


(u, V ) = 


g 1 / 4 \x) 

i 1 / 4 


i — — , 
xy 


1 - — 


xy 


g 


(x 


y ) 1 / 4 


A i-w- 


xy 


(77) 


(78) 


A. 3 Measures 

We recall that the measures are obtained from the direct transitions through 


Res v=u P x \ x (u\v) = 


yx(u) 


(79) 


They have the universal structure 

M x (u) 

Hx(u) = -rexp 

fx(u)fx(-u) 


2K X {u) t ■ M ■ n x (u) — 2k x {uY ■ M ■ k x {u) , (80) 
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with 


M^u) 

M^(u) 


M Fa (u) 


COsh(7T-u) ’ 

7 \g b t A u 

—i - , , 

yjx sinh(7ra) x 2 — g 2 

(-i)Vr(i + f + m)r(i +1 -iu) 

r(a)(U +a W _a ] — g 2 ) a/'(( aJ +a ]) 2 — g 2 )((x^~ a l) 2 — g 2 ) ’ 


for a > 0 . 


(81) 


Moreover, we have that y, F _ a (u) = n Fa (u ) and n^(u) = —/x^(u). Upon analytical continua¬ 
tion to the small fermion sheet, we obtain 


M 0 s (m) 


\Jx 2 — g 2 


(82) 


A.4 Zero momentum limit 


Given a transition P^ Y (u\v) it is immediate to derive its scalings at u = oo, i.e. for a zero 
momentum fermion. This one can be read directly from the function F^ sY (u, v) listed before, 
since the remaining factors in (64) all go to 1 in this limit. We get this way 


P^(u\v) 

g ’ 



U Z/4yl/* ’ 


P^(u\v) 

1 

rs./ - 

\fu ’ 



Pj\i,(u\v) 

~ 1 ( 
g^{ 

y \ 4 / 4 
uJ 


P $\fM V '. 

\ - , 

y/g 

for a < 0 , 

(■ y [+a 

]y[ —a ])l/4 7 

p $\fM v '. 

x ( y [+a 

]y[-a])l/4 

for a > 0 , 

) - 

V9 ’ 

v) and y{ ±a ^ 

= x(v ± 

% |). This information 


(83) 


non-MHV form factors h Y (v) of the excitation Y(v) in the bulk of the paper, with help of 
the asymptotic behaviour of the Jacobian factor 


itself following from 



du 

dPj 


(u) 


r^j 


w 3 / 4 

g 3 / 8 ’ 




ig 1 ' 4 


dpj 

du 


(u) 


r\j 


r 

x cusp 

2 U 2 


(84) 


(85) 
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B The superconformal charge Q and the flux Gold- 
stone fermion 

In section 2.2 we used the realization of the superconformal generator Qa as a zero momen¬ 
tum fermion. The precise relation is (24) and is repeated here for convenience 

s|0> = |p >=. ( 86 ) 
In this appendix we will derive this relation. 

B.l The zero momentum fermion 

The non-trivial part of (86) is the factor dressing the zero momentum fermion state. To 
derive it, we should first fix the normalization of that flux-tube state It is instructive to 
do this in two steps. First we note that we have a well defined flux tube square measure p, 
which allows us to overlap states in the flux Hilbert space. Using the supersymmetry algebra, 
this measure leads to a precise representation of the superconformal generators on the flux, 
that we denote by Qa- There is no reason however for this realization of the supercharge on 
the flux tube to be normalized in the same way as its realization on the generating function 
of amplitudes or equivalently, the super loop. Namely, the two may differ by an overall 
proportionality constant Qa = Co x Qa- 22 We shall now first relate Qa to a zero momentum 
fermion using the measure and the supersymmetry algebra. We will then fix the constant Co 
by demanding that pentagon NHMV amplitude is the same as the MHV one as appears in 
the generating function (14). 

Consider first a delta-function normalized momentum states 

(j)(u)\p(v)) — 2nS(p(u) — piy)) =>- (p = 0\p = 0) = 27t5(0) = VoI(ct) (87) 

where VoI(ct) is the infinite volume of the flux in the coordinate cr conjugate to p. These 
momentum states differ by a simple normalization factor (involving the measure p^,) from 
the rapidity states, which we conventionally normalize as [2] 


IV 27T . . 


, , v, / dv , . , 

\u ) = - — u ) 

Pj{u) 


lp(,l)) _ v 


On the other hand, with respect to this square measure, we have 

2||Qu|0)|| 2 = (0|{Qa,Q a }|0) (89) 

where |0) is the GKP vacuum, normalized so that (0|0) = 1. The commutator is a special 
conformal generator that can be written in terms of the symmetries of the square as 

4 

J]{Qa, Q A } = 2(d r - id 0 ) + C , (90) 

A= 1 

22 See [24] for a very similar relation. 
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where the total hclicity C — 0 in our case and for simplicity, we have summed over the 
R-charged index, see section B.2. The GKP vacuum does not carry U(l) charge while d T 
measures its energy. We conclude that 

||<QU|0)|| 2 = -i?GKp(0|0) = -T cusp Vol(cr) (91) 

where in the last step we used the interpretation of T cusp as the energy density of the flux in 
the a direction [11]. It then follows that 


Q|0) = Um Ip) = ito m> (92) 

where the check over the fermion rapidity (v) indicates that it is on the so-called small 
fermion sheet where the zero momentum point is, see [3]. 

Next, we shall fix the proportionality constant Co- This is done by demanding that the 
pentagon NHMV amplitude is the same as the MHV one. Translated to the POPE notations, 
this condition reads 

1 = P(0|0)=P [ 4 ] (0|0) (93) 


cusp 



2 4 

n 

3 = 1 


lim J-ip- 
dpj, 


^(Vj) X P$4| O (ui,U2,tf3,u 4 |0) x (matrix part) = 


where in the last step we used the large v behaviours quoted in A.4, together with 


g\(u - v) ’ 

that can be read from the expressions in A. 2 , and, finally, the expression for the matrix part 
(which is nontrivial in this case) given by [ 8 ] 

1 

matrix part = —-. 

11 (vi - Vj + i) 
i>j 



We deduce that Co = y/g/Z and hence the relation (86). In the following subsection we 
elaborate on the commutation relation (90). 


B.2 The commutator of superconformal charges 

We shall now derive the relation relation (90) used above. For that aim, it is convenient to 
decompose any twistor in the basis of the square four twistors, (see for example appendix A 
of [ 2 ] for an explicate choice) 

Z = Zfj ^bottom T Zt Z t0 p T Z r bright T Z[ Z} e ft . (94) 
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In this basis, the three symmetries of the square are generated by 
d T z b d Zb ^t^zt i ^rdz r zid Zl and d ^ — ^bd Zb T ^ td Zt z r d Zr %id Z i) • (95) 

The relation (90) is an algebra relation between superconformal generators and therefore 
we can use any representation of the generators to test it. When acting on the generating 
function of helicity amplitudes, the supercharge is represented as Q\ = Z a d v . The operator 
Qa is a specific component of the superconformal generator Q\ which was specified in [9]. 
To translate between (94) and the notations of [9] we may think of the square here as the 
bottom square of the j’th pentagon in the POPE decomposition. Then, equations (10)-(11) 
in [9] for the component of Q\ read 

Qa = d x A oc (Zj_ i A Zj A Zj + 1 ) • Zd v oc z b d v A (96) 

where we used that Zj = Z left , Z 3 _ x = Z right and Z ]+x is a linear combination of Z rig ht and 
Z top . Here, we drop the proportionality factors in (96) as it drops out in the commutator (90). 
We can now use the conjugate operator in this representation to evaluate the commutator 
in (90). We find 23 


Eia.s- 4 } 


yy ; -vV'.'/'v,} = 4 z b d,»+^2 av (a 

A A 

2(d r - idf) + {z b d Zh + z t d Zt + z r d Zr + z t d zi ) + r] A d v A 

A 

2(d T - idf) + C 


Here, the summation over the R-charge index was done for simplicity. Otherwise, on the 
right hand side we would also had an R-charge generator. Alternatively to this derivation, 
(90) can be read from equation (3.9) in [25] by specifying to the corresponding component. 


C A second example: V\ o V 2 0 Vu 

In this second example we analyse another heptagon component where the middle pentagon 
is charged, namely VioV 2 °Vsa- We consider again the contribution from multiparticle states 
in both squares, this time a term proportional to 

e —2n—Sifc/2 x e -2r 2 +i02 _ (93) 

To End which processes contribute we first note that we need twist 2 states for both squares 
with helicities -3/2 and +1. Then we follow the same logic as in section 3.2 and see the 
implications of R-charge conservation. The first pentagon carries one unit of R— charge 
which enforces the state propagating in the first middle square to be in the anti-fundamental 
(4) representation of SU( 4). The second pentagon also possesses one unit of R— charge. This 

23 Note that the commutation relation (97) is independent of the measure one uses to realise it and thus 
{Qu,Q a } = {Qu,Q a }. 
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implies that the state in the second middle square should be in one of the two representations 
entering in the decomposition of 4 <g> 4 = 6 © 10. The state is finally projected onto the 
vector representation (6) by the last pentagon that carries two units of R— charge. In the 
first square, the only option is the state formed by 0R_i. In the second one we could either 
have 4>Fi or ' 0 - 0 ; as the reader might guess, it is the second state which contributes at tree 
level given the small fermions behaviour. All we are left to do is study the process 

vacuum —> 0 F _i —> 00 —>• vacuum. (99) 

We see that the R-charge index of the antifermion in the first middle square is fixed by the 
R.-charge index of x hi the first pentagon. However, this is not the case for the two fermions 
in the second square. In other words the matrix part is not trivial in this case, but one can 
compute it in a straightforward way following the logic in [3] where the contribution of (00) 


was studied and the S-matrix of the fermions in the vector (6) representation [5]. 
integrand is then given by the product of the following parts 

The full 

dynamical part = 

VF-Au) /ty(w 2 ) Mv) 

(100) 


w P 4<\f 1 (w2\u)P^ Fi (v\u)P^(w 2 \ w i)P^(v\w 1 ) 
P Fl ^(u\w 1 )P^ F _ 1 (w 1 \u)P^(w 2 \v)P^(v\w 2 ) 1 

(101) 

form factors part = 

\h^{wi)h^(wi)h F _ 1 (u)h^(w2)h^{v){h i ,{w 2 )) 2 {h^(v))‘ 2 

9 3 

matrix part = 

2 

(103) 

(v - w 2 ) 2 + 1 ’ 


Given that the fermions can be small or large, one can split the process into the different 
contributions 


F f _]_—F F_ 1 ^q/, s i/j s T 2 x VV^ S > (105) 

where the factor of 2 accounts for the two equivalent ways of choosing which of the two 
fermions is small. At tree level only the last term contributes. One could expect that also 
the second to last term with two small fermions in the last square would also appear at tree 
level, but as explained in [3] this contribution vanishes 24 . Therefore we consider only the 
last term in (104) and the transition reads 

f du f dv f dw 2 f dwi . 

** = 2 J J 2* J ./ ^7 (Vnamical part) 

R R—iO Csmall Csmall 

x(form factors part) x (matrix part), (106) 

where we use —ie prescription for fermions in the top square and the opposite prescription for 
fermions in the bottom square. (Accordingly we use the lower half plane (half-moon) contour 

24 This is because after integrating out the first small fermion there are no more singularities enclosed by 
the second integration contour. In colloquial words, one should have at least one large excitation to attach 
the fist small fermion (e.g. Fibs0s would have a non vanishing contribution). 
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Csmaii for the small fermion at the bottom and its conjugate C sma u for the small fermion at 
the top.) We can now integrate out the two small fermions by picking the poles w-i — v + i 
coming from the pentagon transition and W\ = u — i/2 from the matrix part. Finally 

we integrate over u and v and find that the tree-level prediction is given by 


hk/ ; J7 _ j 


g—2 ti —2t 2—3i0i /2+i(j>2 


( e 2 <ri _|_ 2 ) 2 (e 2a2 + 1) 2 (e 2<Tl + e 2<T2 + e 2ai+2a2 ) 3 
x [e 7ai + 3e 5ffl+2<T2 + 5e 7ai+2a2 + 3e 3ai+4tT2 + I2e 5ai+4a2 
+ 10e 7ai+4<72 + e ffl+6fT2 + 5e 3ai+6a2 + 10e 5ffl+6ff2 + 6e 7ai+6a2 ) . (107) 


To check this expression against data, we first relate TioT 2 o V 34 to a linear combination of 
super amplitude //-components [9] 


d d 


dx\ dx.2 \dx3 


d 


W = 


-l)i ((5 ) e 


d 


(1) 2 (2) 2 (3) 2 dr }_i 


d 


d 


{w '~ 1) 9 n X + {w ' 0) WJ w 


d 


( (iIm[ 5 (3 3 )1 ^ 2 ’ 3 ’ _1) W(_1 " 1,5 ’ 5) + 2 ’ 3 ’ °) W ( ' 1,0 ’ 5 ’ 5) ) (108) 


Then we extract the amplitude components in Mathematica and evaluate the full expression 
with the heptagon twistors. After expanding at large T\ and r 2 and picking up the term 
proportional to (98) we find again a perfect match with the OPE result (107). 
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